Estabilidad y dinámica de modos ópticos en estructuras no lineales, helicoidales y PT-simétricas

Dr. Claudia Castro Castro
2 de Diciembre del 2020

Seminario “Dr. Alberto Rubio” Facultad de Ciencias UABC

Índice

  • Introducción
    • Definición de onda
    • Antecedentes
    • Fibras ópticas
    • Fenómenos ópticos lineales vs no lineales
    • Ecuaciones de Maxwell
    • ¿Qué significa PT-simétrico?
    • ¿Qué significa torcido?
  • Ecuación de Schrödinger
  • Dinámica y estabilidad
  • Conclusiones

Definición intuitiva de onda

    Cualquier señal reconocible que se transfiere de una parte del medio a otra con una velocidad de propagación reconocible

Antecedentes

1790s: “Telégrafo óptico” inventado por Claude Chappe Telegrafo optico y sus simbolos

Courtesía de Wikipedia Commons

Luz guida en chorros de agua

  • 1840s: Daniel Collodon y Jacques Babinet demostraron que la luz podía guiarse a través de chorros de agua para exhibiciones de fuentes.

  • Cuando los rayos de luz en el agua golpean el borde del chorro con un ángulo, la reflexión interna total los atrapa en el líquido.

  • Refracción en la superficie depende de la diferencia en el índice de refracción; cuanto mayor es la diferencia, más refracción

Reflexión interna total (TIR)

  • Es un efecto secundario de la refracción, la curvatura de la luz que pasa de un material transparente a otro.
Refraction de vidrio a aire

Courtesía de Hecht, J. City of Light: The Story of Fiber Optics, 1999


  • Índice de refracción \[ n=\frac{\text{velocidad de la luz en el vacio}}{\text{velocidad en el material}} \]

  • Ley de Snell \[ n_1 \sin{\theta_1} = n_2 \sin{\theta_2} \]

  • Resposable por
Lente de faro Fresnel

Flor de loto de cristal

Total Internal Reflections in Diamond

Courtesía de Wikimedia Commons

Reflexión interna total (TIR)

Luz guida en chorros de agua

Courtesía de Hecht, J. City of Light: The Story of Fiber Optics, 1999

  • La luz que golpea la superficie casi directamente se refracta en el aire
  • Hay un ángulo crítico donde no puede emerger al aire, se mide a partir de una línea llamada normal que es perpendicular a la superficie.
  • Si \( \theta_1 > \sin^{-1}\left(\frac{n_2}{n_1}\right) \) \( \Rightarrow \sin \theta_2 > 1 \;\;(!) \)
  • La luz fuera del ángulo crítico se refleja de nuevo en el cristal.
  • ángulo crítico \( \theta_c = \sin^{-1}\left(\frac{n_2}{n_1}\right) \) existe solo cuando \( n_2 < n_1 \), entonces TIR ocurre solo dentro del medio con índice más alto
Demonstration of Total-Internal-Reflection(TIR) in a wine glass

Courtesía de Keerthi, CC BY 4.0 via Wikimedia Commons

Inicios de las fibras ópticas

  • Generaciones de artesanos, inventores, ingenieros y científicos desarrollaron la tecnología para hacer vidrio
    • Es maleable cuando está caliente, así como la transparencia del sólido, lo que lo ha convertido en un material atractivo.
    • El ingrediente más importante del vidrio común es el dióxido de silicio, un mineral duradero conocido como Sílice/Silica.
  • 1930s: Heinrich Lamm (médico) primera persona conocida que ha demostrado la transmisión de imágenes a través de un arreglo de fibras ópticas
    • Sin embargo, las fibras descubiertas transmitían imágenes de manera deficiente.
  • 1950s: Van Heel y O'Brien recubren una fibra descubierta con un material de índice menor mantendría la TIR mientras protege la superficie óptica
    • Fibras revestidas de vidrio tenían una atenuación de aproximadamente un decibel por metro (dB/m), lo que está bien para imágenes médicas, pero demasiado alto para comunicaciones.

Características de las fibras ópticas

Esquema de fibra optica multi-nucleo

Courtesía de depositphotos.com

Esquema de un arreglo 2D

Pertsch et al. Nonlinearity and disorder in fiber arrays. Physical Review letters, 2004

Esquema de un arreglo 2D

Joannopoulos et al. Photonic crystals: molding the flow of light. Princeton University Press, 2011.

Características de las fibras ópticas

Esquema de fibra optica dielectric y el perfil de indices de refraccion

Agrawal G., Nonlinear fiber optics, Academic Press, 2013

Dispersión

Gota que cae en agua

Cortesia de Pixabay

  • Como cuando una piedra que cae en un estanque inmóvil produce ondas en la superficie del agua, eventualmente desaparecen y se desvanecen
  • Un rayo de luz puede extenderse cuando viaja a través de diferentes medios.
  • La dispersión temporal y la difracción espacial se deben a la dependencia del índice de refracción de la frecuencia y la longitud de onda
  • Ambos causan un ensanchamiento espacial y temporal de la luz.

Ecuaciones de Maxwell

La luz es una onda electro-magnética que consiste de un campo eléctrico y un campo magnético oscilando a una tasa muy alta (\( 10^{14} \) Hz) viajando en el espacio con una dirección perpendicular a ambos campos vectoriales

  • \[ \begin{align} \nabla\times\mathbf{E}=&-\frac{\partial\mathbf{B}}{\partial t}\\ \nabla\times\mathbf{H}=&\mathbf{J}+\frac{\partial\mathbf{D}}{\partial t}\\ \nabla\cdot\mathbf{D}=&\rho_{v}\\ \nabla\cdot\mathbf{B}=&0 \end{align} \]

  • \( \mathbf{E} \) campo eléctrico [\( V/m \)]
    \( \mathbf{B} \) densidad de flujo magnético [\( T \)]
    \( \mathbf{H} \) campo magnético [\( A/m \)]
    \( \mathbf{D} \) densidad de flujo eléctrico [\( C/m^2 \)]
    \( \mathbf{J} \) densidad de corriente eléctrica [\( A/m^2 \)]

  • \( \mathbf{D}=\varepsilon \mathbf{E} \)
    \( \mathbf{B}=\mu \mathbf{H} \)
    \( \mathbf{J}=\sigma \mathbf{E} \)

    \( \varepsilon= \varepsilon_0 \varepsilon_r \) dielectric permittivity
    \( \mu= \mu_0 \mu_r \) permeability
    \( \sigma \) electric conductivity

    Ecuación de onda

  • El campo vectorial eléctrico satisface la ecuación de onda \[ \nabla^2 \mathbf{E}(\mathbf{r},t)-\frac{1}{c^2}\frac{\partial^2 \mathbf{E}(\mathbf{r},t) }{\partial t^2}=0 \]
  • T-simetría: ecuación es invariante w. r. t. tiempo \( t\rightarrow -t \)

Óptica lineal vs no lineal

  • Lineal : las ondas electromagnéticas inducen una separación de las cargas en el material, es decir, una polarización \( P_L \), la cual es directamente proporcional al campo eléctrico \[ P = \varepsilon_0 \chi^{(1)} E \] donde \( \varepsilon_0 \) es la permitividad del vacío y \( \chi^{(1)} \) es la constante de susceptibilidad.
  • Interacción de la luz con la materia no modifica las propiedades de onda
  • Ejemplos de fenómenos ópticos lineales:
Reflejo de un gato en un espejo

Reflexión

Ejemplos de fenómenos ópticos lineales

Refracción

Lentes photocromaticos

Photocromia

  • \( \diamond \) Linealidad es una suposición que solo es válida para bajas intensidades

  • \( \diamond \) Casi todos los materiales tienen algunos efectos no lineales si la fuente de luz solo es lo suficientemente potente altas intensidades

Óptica no lineal

  • No lineal: describe el comportamiento de la luz en medios en cual el componente dieléctrico de la polarización responde a la forma no lineal del campo eléctrico de la luz \( E \) \[ P=\varepsilon_0 \chi^{(1)}E + \varepsilon_0 \chi^{(2)}E^2 + \varepsilon_0 \chi^{(3)}E^3+\dots \]
  • 2nd order susceptibility \[ \chi^{(2)} \approx 1.94\times10^{-12}\;m/V \]
  • 3rd order susceptibility\[ \chi^{(3)} \approx 3.78\times10^{-24}\;m^2/V^2 \]
  • Para fibras de sílice \( SiO_2 \) el effecto de segundo orden es despreciable
  • \( \chi^{(3)} \) está asociada al effecto Kerr
  • Polarizacion total: lineal + no lineal \[ P(\omega) = \varepsilon_0 \chi^{(1)}E(\omega)+3\varepsilon_0\chi^{(3)}|E(\omega)|^2E(\omega) \]
  • Effective susceptibility \[ \chi_{eff} = \chi^{(1)}+3\chi^{(3)}|E(\omega)|^2 \]
  • Está ligada al índice de refracción \[ \bar{n} = 1+\chi^{(3)}=\bar{n}_0 +\bar{n}_2\;I \]

Óptica no lineal: effecto Kerr

  • Descubierta por John Kerr en 1875
  • Describe situaciones donde el índice de refracción depende de el campo eléctrico como \[ \bar{n}(\omega, |E|^2)=\bar{n}_0(\omega) + \bar{n}_2(\omega)|E|^2 \]
  • Para sílice es \( 1.3\times10^{-22}\; m^2/V^2 \)

Derivación de la ecuación de Schrödinger

  • En medio Kerr el índice de refracción depende de la intensidad del campo eléctrico \( I(t) \) \[ \bar{n}(t)=\bar{n}_0 + \bar{n}_2 I(t), \;\;\;I(t)=2\bar{n}_0\varepsilon_0 c |A(z,t)|^2 \]
  • y \[ E(z,t)=A(z,t)e^{i(\omega_0 t-\beta _0 z)} \]
  • Aplicando transformada de Fourier del campo óptico \[ E(z,t)=\int_{-\infty}^{\infty} \tilde{E}(z,\omega)e^{i(\omega t - \beta z)} \;d\omega \]
  • Considerar la expansion de Taylor de la constante de propagación \[ \beta(\omega) = \beta_0 + \beta _1 (\omega - \omega_0 )+\frac{1}{2}\beta_2 (\omega - \omega_0)^2 + \Delta \beta_{NL} \]
  • Sustituir \[ \begin{align} E(z,t)=&e^{-i\beta_0 z}\int_{-\infty}^{\infty} \tilde{E}(z,\omega)e^{i\omega t - i\beta_1 z \Delta \omega - i\frac{1}{2}\beta_2 z \Delta \omega^2-iz\Delta\beta_{NL} } \;d(\Delta\omega)\\ =&e^{i(\omega_0t-\beta_0 z)}\int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)e^{it\Delta\omega-i\beta_1 z \Delta \omega - i\frac{1}{2}\beta_2 z \Delta \omega^2-iz\Delta\beta_{NL}} \;d(\Delta\omega)\\ \equiv&e^{i(\omega_0 t-\beta _0 z)} A(z,t) \end{align} \]

Derivación de la ecuación de Schrödinger

  • Hemos definido \[ A(z,t) = \int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)e^{it\Delta\omega-i\beta_1 z \Delta \omega - i\frac{1}{2}\beta_2 z \Delta \omega^2-iz\Delta\beta_{NL}} \;d(\Delta\omega) =\int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)e^{ig(z,t)} \;d(\Delta\omega) \]
  • Derivar \[ \begin{align} \frac{ \partial A(z,t)}{\partial t} =&\int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)i\Delta \omega e^{ig(z,t)} \;d(\Delta\omega) \\ \frac{ \partial^2 A(z,t)}{\partial t^2} =&\int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)(i\Delta \omega)^2 e^{ig(z,t)} \;d(\Delta\omega) \\ \frac{ \partial A(z,t)}{\partial z} =&\int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega)(-i\beta_1 \Delta \omega - i\frac{1}{2}\beta_2 \Delta \omega^2-i\Delta\beta_{NL}) e^{ig(z,t)} \;d(\Delta\omega) \end{align} \]
  • Combinar \[ \frac{ \partial A(z,t)}{\partial z} + \beta_1 \frac{ \partial A(z,t)}{\partial t} -i\frac{1}{2}\beta_2 \frac{ \partial^2 A(z,t)}{\partial t^2} = \int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega) \left[ ... \right]e^{ig(z,t)} \;d(\Delta\omega) \]
  • La expresion en \( [...] \) es \( -i\bar{n}_2k_0 I \)
  • Entonces \[ \frac{ \partial A(z,t)}{\partial z} + \beta_1 \frac{ \partial A(z,t)}{\partial t} -i\frac{1}{2}\beta_2 \frac{ \partial^2 A(z,t)}{\partial t^2} = -i\bar{n}_2k_0I \int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega) e^{ig(z,t)} \;d(\Delta\omega) \]

Derivación de la ecuación de Schrödinger

  • Entonces \[ \frac{ \partial A(z,t)}{\partial z} + \beta_1 \frac{ \partial A(z,t)}{\partial t} -i\frac{1}{2}\beta_2 \frac{ \partial^2 A(z,t)}{\partial t^2} = -i\bar{n}_2k_0 I \int_{-\infty}^{\infty} \tilde{E}(z,\omega_0 +\Delta \omega) e^{ig(z,t)} \;d(\Delta\omega) \]
  • \[ ``\qquad" \qquad = -i\bar{n}_2k_0 I A(z,t) \]
  • Por lo tanto
  • \[ \frac{ \partial A(z,t)}{\partial z} + \beta_1 \frac{ \partial A(z,t)}{\partial t} -i\frac{1}{2}\beta_2 \frac{ \partial^2 A(z,t)}{\partial t^2} = i\gamma |A(z,t)|^2A(z,t)-\frac{\alpha}{2}A(z,t) \]

Ecuación discreta no lineal de Schrödinger (DNLS)

  • Considere la dinámica de propagación en un arreglo discreto de fibras ópticas descrita por la ecuación discreta de Schrödinger no lineal (DNLS) \[ \begin{equation}\label{eq:DNLSE} i \frac{dc_n}{dz} = \epsilon_n c_n - k(c_{n+1}+c_{n-1}) + \sigma|c_n|^2c_n, \end{equation} \]
  • \( c_n \): amplitudes complejas, en el sitio \( n \), que dependen de la dirección de la variable de propagación \( z\in\mathbb{R} \)
    \( k \): fuerza uniforme de acoplamiento con el vecino más cercano
    \( \epsilon_n \): perfil de índice de refracción in-situ
    \( \sigma \): fuerza de no linealidad
  • \( \diamond \) El hamiltoniano que da lugar a las ecuaciones de movimiento viene dado por la ecuación \[ \begin{equation}\label{eq:hamiltonian} H_D = \sum_n \epsilon_n|c_n|^2 + \frac{\sigma}{2}|c_n|^4 - k\left(c_{n+1}c_n^*+c^*_{n+1}c_n \right). \end{equation} \]
  • Configuración ideal de guías de ondas ópticas idénticas, donde \( \epsilon_n \) es una constante fija \( \epsilon \).
Surface with positive orientation

Representación esquemática de un arreglo de guías de ondas ópticas planar. \( r_n \) es el radio de cada núcleo en particular

Parity-Time (PT) simetría

  • Quantum Mechanics Optics
    Schrödinger equation Paraxial equation
    \( i\hslash\frac{\partial\Psi}{dt}+\frac{h^{2}}{2m}\frac{\partial^{2}\Psi}{\partial x^{2}}-V\left(x\right)\Psi=0 \) \( i\lambda\frac{\partial E}{\partial z}+\frac{\lambda^{2}}{2}\frac{\partial^{2}E}{\partial x^{2}}+n\left(x\right)E=0 \)
    time propagation distance
    \( t \) \( z \)
    Plank’s constant Wavelength
    \( \hslash \) \( \lambda=\frac{1}{k} \)
    Probability amplitude Electric field envelope
    \( \Psi\left(x,t\right) \) \( E\left(x,z\right) \)
    Complex potential Complex refraction
    \( V\left(x\right)=V_{R}\left(x\right)+iV_{I}\left(x\right) \) \( n\left(x\right)=n_{R}\left(x\right)+in_{I}\left(x\right) \)

Parity-Time (PT) simetría

  • Bender consideró si un sistema hamiltoniano mecánico cuántico con un potencial complejo puede tener un espectro real.
  • P: \( \hat{x}\rightarrow -\hat{x}\; \) y \( \;\hat{p}\rightarrow -\hat{p} \)
    T: \( \hat{x}\rightarrow \hat{x}\; \), \( \;\hat{p}\rightarrow -\hat{p}\; \), y \( \;i\rightarrow -i \)

    \( \hat{p} \): momentum operator.
    \( \hat{x} \): space operator.
  • Condición de PT-simetría \[ V(x) = V^*(-x)\qquad n(x)= n^*(-x) \]

    Esto se traduce en el potencial complejo cuya parte real es una función par mientras que la parte imaginaria es impar.
  • \( \circ \) Hamiltoniano PT-simétrico
    • Puede tener un eigen-espectro real
    • Eigen-estados PT-simétricos puede que occurran para cierto régimen paramétrico

Parity-Time (PT) simetría

Sistemas ópticos acoplados convencionales y PT-simétricos.

Rüter et al., Nature Physics, 2010

Propagación del haz en dos guías de ondas PT-simétricas no lineales. Metasuperficie

Ramezani et al. PR A 2010

Dispersor recubierto con una metasuperficie ultrafina con pérdida y ganancia equilibradas. Metasuperficie

Sounas et al., Phys. Rev. Applied 2015

Fibras multi-núcleo no lineal PT-simétricas

  • Considere la dinámica de propagación del haz en una fibra discreta de múltiples núcleos de \( N \) sitios dispuestos igualmente espaciados en un anillo de radio \( R_0 \) donde cada núcleo tiene un radio \( r_0 \), descrito por las ecuaciones de modos acoplados Dinámica de propagación en una matriz discreta de fibras ópticas descrita por la ecuación discreta de Schrodinger no lineal (DNLS) \[ i\frac{{dc_{n}}}{dz}=k\left(c_{n+1}+c_{n-1}\right)+i\gamma_{n}c_{n}+\sigma|c_{n}|^{2}c_{n} \]
  • \( c_n \): representa amplitudes complejas, en el sitio \( n \), que dependen de la dirección de la variable de propagación \( z\in\mathbb{R} \)
    \( k \): fuerza uniforme del acoplamiento con el vecino más cercano
    \( \gamma_n \):tasa de ganancia óptica (\( {\color{blue} {\gamma_n>0}} \)) o pérdida (\( {\color{gray} {\gamma_n<0}} \))
Esquema de una fibra optica torcida

Esquema de fibra óptica multi-core PT-simétrica

Fibras multi-núcleo no lineal helicoidal PT-simétricas

  • Considere la dinámica de propagación del haz en una fibra discreta de múltiples núcleos de \( N \) sitios dispuestos igualmente espaciados en un anillo de radio \( R_0 \) donde cada núcleo tiene un radio \( r_0 \), descrito por las ecuaciones de modos acoplados Dinámica de propagación en una matriz discreta de fibras ópticas descrita por la ecuación discreta de Schrodinger no lineal (DNLS) \[ i\frac{{dc_{n}}}{dz}=k\left({\color{orange} {e^{-i\phi}}}c_{n+1}+{\color{orange} {e^{i\phi}}}c_{n-1}\right)+i\gamma_{n}c_{n}+\sigma|c_{n}|^{2}c_{n} \]
  • \( c_n \): representa amplitudes complejas, en el sitio \( n \), que dependen de la dirección de la variable de propagación \( z\in\mathbb{R} \)
    \( k \): fuerza uniforme del acoplamiento con el vecino más cercano
    \( \gamma_n \): tasa de ganancia óptica (\( {\color{blue} {\gamma_n>0}} \)) o pérdida (\( {\color{gray} {\gamma_n<0}} \))
    \( \phi \): Peierls phase
Esquema de una fibra optica torcida

Esquema de fibra óptica helicoidal PT-simétrica

Ruptura de PT-simetría: caso lineal

\[ i\frac{{dc_{n}}}{dz}=k\left(^{-i\phi}c_{n+1}+^{i\phi}c_{n-1}\right)+i\gamma_{n}c_{n} \]
  • PT-simetría alternante \( \gamma_n=(-1)^{n+1}\gamma \), \( \gamma=2.04 \)
    • Cambio de signo \( \Leftrightarrow \) ganancia/perdida
  • Condiciones de frontera periódicas \( c_{n + N} = c_n \)
  • Otros valores usados \( \sigma=0 \), \( n_s=1.552 \), \( R_0=8\mu m \), \( \lambda=980nm \), \( r_0=0.91R_0 \)
Observe the PT-symmetry breaking induced by the twist of the fiber. Esquema de una fibra óptica torcida

Longhi S., PT phase control in circular multi-core fibers. Optics letters, 2016.

  • Cuando se tuerce por un ángulo de \( \phi =\pi/6 \) se volve inestable :(

No lineal pero conservado

\[ i\frac{{dc_{n}}}{dz}=k\left(e^{-i\phi}c_{n+1}+e^{i\phi}c_{n-1}\right)+{\color{gray} {i\gamma_{n}c_{n}}}+\sigma|c_{n}|^{2}c_{n} \] con \( \gamma_n =0 \)

  • If we make the transformation \( c_n = a_n(z)e^{i\varphi_ n} \), \( \varphi_ {n+1} −\varphi_n= \phi \)
  • The equations for \( a_n \) are those of the classical uniform array, but now the periodic boundary condition is modified to \( a_{n+N} = a_n e^{ iN\varphi} \).
  • This transformation also shows in the Hamiltonian \[ H=\sum_{n=1}^N k \left( c_{n + 1} c_n^*e^{-i \phi} + c_n c^*_{n+1} e^{i \phi}\right)+\frac{\sigma}{2}|c_n|^4 = \sum_{n=1}^N k \left( a_{n + 1} a_n^* + a_n a^*_{n+1} \right)+\frac{\sigma}{2}|a_n|^4. \]

Conclusiones

  • Exploramos cómo la dinámica de la fibra no lineal con seis núcleos se ve afectada por un giro inducido en el escenario cuando la ganancia/ pérdida no está presente en el modelo y cuando hay un perfil de ganancia/ pérdida alternante.
  • Los resultados numéricos destacan el potencial de inducir una torsión de fibra para controlar la dinámica de la luz en fibras de múltiples núcleos no lineales y sugieren un escenario rico para una mayor exploración del espacio de parámetros.

  • ¡Gracias!